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Tempered distributions: does universal 
tempering procedure exist? 
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Abstract 

Since the turn of the century, there has been increased interest in 
the application of heavy-tailed distributions, particularly stable dis¬ 
tributions, to problems in physics and finance. Although, the tails of 
stable distributions provide a better fit to real-world data, they are 
too fat to describe empirical distributions. To remedy this drawback 
of stable distributions, so-called tempered variants of stable distribu¬ 
tions have been proposed. In this paper, we argue that the tempering 
should be connected to the model leading to heavy tailed distribution 
and propose several tempering procedures in connection with corre¬ 
sponding models 

Keywords: stable distributions; tempered stable distributions; LePage 
series; mixtures of normal distributions; selling-short-strategy 

1 Introduction 

Stable distributed are widely used in many branches of the sciences and 
economics because of the overwhelming empirical evidence that real-world 
distributions exhibit fatter tails than that of the Gaussian distribution. Ex¬ 
amples in financial economics include the works of Mandelbrot [1], Fama 
[2], Samuelson [3], Embrechts et al. [4], Rachev and Mittnik [5j. However, 
it is clear, that all amount of money in the world is finite. This implies 
paradoxically that heavy tailed distributions may provide a good model for 
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large, but not too large money sums. Accordingly, one may expect, that the 
tails of stable distributions are too heavy for consideration of large money 
amounts. Connection of this idea with applicability of generalized limit the¬ 
orem was proposed in [am. Despite the heavier tails than the Gaussian 
distribution, stable distributions have been found too heavy than found in 
real-world distributions (see, for example, H31). The most popular approach 
for correcting for heavier tails than warranted is the so-called tempering 
procedure that involves modifying the characteristic function of stable dis¬ 
tribution by changing of Levy measure in Levy-Khinchine representation of 
stable distributions. The characteristic function ficrs f° r a tempered stable 
distribution is given by 

4>cts(u) = exp^w/x+C'iT(—a)((A + —m)“—A“)+C 2 T(—a)((A_+m)“—A“)^, 
for some p G IR. Moreover, ficrs can be extended to the region 

{z;6C: fm(^) <G (—A_,A + )}. 

Despite the popularity of this approach, we demonstrate that this procedure 
is not universal and any modification of the stable distribution should be 
made on an individual basis. We also are that for non-stable distributions 
with heavy tails, the tempering procedure makes no sense for such distribu¬ 
tions. 


2 Examples of models leading to heavy tailed 
distributions 

Here we give examples of toy-models (i.e. small models ignoring many details 
of the process studied, but taking into account one of the most essential part 
of it) for probability distributions with heavy tails. 

Example 2.1. The First Passage Time Distribution. Here we are 
interested in the probability that a particle reaches a given point c at a time 
t. It is well-known that first-passage time for a Brownian particle follows a 
Levy distribution (see W- This distribution is |- stable. It is concentrated 
on the positive semi-axis and has characteristic function 

fit ) = exp{ —\/— 2 bit}, t G IR 1 , b > 0. (2.1) 

This distribution has a heavy tail of order 1/2. 
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Example 2.2. Discrete variant of Example \2.1[ Now consider the 
random walk on real line with equal probabilities of moving to the right or 
to the left. We are interested in the distribution of the time of the first 
passage through 1. It is known that the probability generating function of this 
distribution is 

m (2.2) 

Clearly, the corresponding distribution has a heavy tail of order 1/2. 

Example 2.3. Sub-Gaussian distributions. Let X be a Gaussian ran¬ 
dom variable with zero mean. Suppose that A is a-stable positive random 
variable (it is necessary to have a G (0,1), and skewness parameter (5 = 1). 
If X and A are independent, then the product A 1 / 2 X has a symmetric 2a- 
stable distribution. This is the so-called sub-Gaussian distribution. Clearly, 
the sub-Gaussian distribution can be represented as a product of two stable 
distributed random variables. Such representation is not unique. 


Example 2.4. LePage series. The LePage series allows to obtain many 
examples of toy-models leading to heavy tailed distributions (see, for example, 
m > • Suppose, we have a Poisson process, and its arrival times are Tg, j = 
1,2,.... Let {Xl, X- 2 ,...} be a sequence of independent identically distributed 
(i.i.d.) random variables and suppose that this sequence is independent with 
Poisson process mentioned above. Consider the following series 


OO 


£ 



(2.3) 


where a is a real number between 0 and 2 (strictly), and random variable Xi 
has an absolute moment of order greater that a. If the series 12.31) converges, 
then its sum follows a strictly stable distribution with index of stability equals 
to a. 


The sum of the series (12.31) allows for the following three interesting in¬ 
terpretations. 


1 . 


Consider electrical charges Xj, (j = 1,2,...) on a straight line, ar¬ 
ranged with accordance to the distances X, from the origin, where is a 
fixed charge Q. According to Coulomb’s law, the force acting on the 
charge Q is 


OO 


i=i 
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Under condition that all charges Xj are chosen as i.i.d. random vari¬ 
ables, we obtain that the random force Y has a strictly stable distri¬ 
bution with the parameter a = 1/2. 

2. Similar to the first interpretation suppose that instead of charges we 
have masses, and use Newton law instead that of Coulomb’s law. Again, 
the force acting on the mass Q will have a strictly stable distribution 
with stability parameter a = 1/2. In this situation, all Xj are positive, 
and, therefore, random variable Y is positive, too. Now we may con¬ 
clude, that Y has Levy distribution. Note that we may consider the 
Poisson held on Euclidean plane or in Euclidean space. The distances 
T 3 from a particle to the origin can be considered as the arrival times 
of a Poisson process. Therefore, we will have the same interpretations 
for the mass (or charges) distributed according to the Poisson held. 

3. Consider Euclidean plain IR 2 , and suppose that there is a mobile phone 
base station at the origin. Suppose further that mobile phones are 
distributed on the plane according to a Poisson held. The jth phone 
sends signal Xj to the base station. The signal decays inversely to some 
exponent of the distance between the phone and base station (i.e., the 
signal, coming to the station from a phone is Xj/T^ 01 ). From technical 
experiments it is known, that 1/a is a number around 2.6. Now we can 
say that the summand signal coming to the base station has strictly 
stable distribution. 

Example 2.5. Pareto distribution. Let S = {AR,A" 2 ,...} be a sequence 
of i.i.d. positive random variables. Suppose also that {u p , p e (0,1)} is a 
family of random variables, independent of the sequence S and having geo¬ 
metric distribution with parameter p: lP{z/ p = k} = p( 1 — p) k ~ x , k = 1, 2,.... 
Consider the normalized product of random number of the elements of S, that 


is, 


z P = U x 7- 


(2.4) 


It is known (see m) that the limit distribution of Z p as p —> oo is a Pareto 
distribution with shape parameter a = 1 / 7 , where 7 = IElogAR, assuming 
this parameter is positive, that is, 



forx > 1 , 
forx < 1 , 


(2.5) 
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in the case when 7 > 0. This particular case is interesting for us because it 
leads to heavy tailed distribution. 

Also note that if random variable X\ has a Pareto distribution with pa¬ 
rameter a > 0, then 


Up 



3 = 1 


where X 1 ,X 2 ,... are i.i.d. random variables, and {v p , p e (0,1 )} is a family 
described above. The sign = is used to denote the equality in distribution. 

Note that the distribution of the value of Z p represents the distribution 
of capital, obtained in random market over v p units of time. The mean value 
of v p is large in this situation (or, equivalently, the parameter p is small). 
On economical interpretation of Z p and its limit distribution the reader is 
addressed to HI] 

Example 12. 5 1 shows that the reason for the distribution having heavy tails 
is the randomness of the number of multipliers in the product (12.4p . This 
number can be arbitrarily large with a positive probability, while the multipli¬ 
ers may have thin tails. However, note that we do not have a sum of a random 
number of random variables, but a product of them. In the case of sums it 
is impossible to obtain heavy tailed limit distributions using summands with 
thin tails. Accordingly the limit theorem on convergence to stable distribu¬ 
tion cannot be considered as an explanation for the appearance of heavy tailed 
distributions in applied problems. The central limit theorem shows, we need 
to have the summands with heavy tails for sum to be convergent to stable 
distribution. 

Example 2.6. Selling-short-strategy in financial market. Suppose 
that an investor follows a selling-short-strategy. This means that the investor 
has borrowed some assets and sells it at Price of P* with the expectation that 
the price will decline below. Should the price decline, the investor will buy 
the assets at the lower price and returns the asset to the original owner. 
Let us study this strategy in more detail for the discrete case. Suppose that 
{Pj, j = 1 . 2 ,...} is a sequence of prices for an asset at times 1,2,.... In 
our toy-model we consider Pi, P 2 ,... as a sequence of i.i.d. random variables 
(which is, of course, a simplification). The short-seller borrows the asset, 
and will sell it at time 1 , 2 ,... prior to the first appearance of the event 
{P < P*}. The moment of this appearance will be denoted by v + 1. It is 
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clear that the time elapsed before the onset of the event may be considered a 
random variable. A first approximation for describing this random variable 


is to consider it as a number of experiments in a Bernoulli scheme until the 
first appearance of a success. The probability of success is p = P{P < P*}. 
In other words, the random number of experiments has geometric distribution 
(starting at 1) with parameter p 6 (0,1). Now we see that the sum obtained 
by the short-seller in this strategy is 



(2.6) 


Here Pj is the price at the time j, X 3 is the quantity of assets sold by the short- 
seller at this time, and v p is a random variable with geometric distribution 
P{i/ p = k} = p( 1 — p) fc_1 , k = 1,2,.... For simplicity, we suppose that u p is 
independent on both sequences {Pj, j = 1,2,...} and {X,-, j = 1,2 ,..and 
last two sequences are independent on each other as well. Let us make some 
suppositions on the character of the sequence {Xj, j = 1,2,...}. Random 
variable Xj is an integer equals to the number of assets sold at moment 
j. It is more or less natural to consider Xj as a number of successes in 
Bernoulli scheme, but with the probability of fail depending on the number of 
the experiment. Namely, we consider the process of selling at the time j in 
the following way. It is a sell of the first asset, after that a sell of the second 
asset, and so on, until nobody will buy assets. However, the probability to 
sell second asset is supposed to be higher, than for the first asset; probability 
to sell third asset is higher than for the second one; and so on. It can be 
connected to the fact that most gamblers (according to some sources, up to 
80%) prefer strategy “to buy”. For the first approximation, we suppose, that 
the probability to fail while selling kth asset is 7 /k, where 7 is the probability 
to fail while selling the first asset. Under these suppositions, we see that all 
random variables Xj, j = 1,2,... are i.i.d. with the distribution of the form 


k -1 



It is clear that probability generating function of Xj is 

Q(z) = 1 - (1 - z)\ 


(2.7) 
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Let Lp(s) be Laplace transformation of random variable P\. Then Laplace 
transform of the product PjXj is 


oo oo k —1 

Lp X (s) = Eexp {—sPjXj } = L P (sk)W{X j — k} — L P (s£;)^ JJ(i-y). 

k =1 k =1 i =1 

Now it is easy to calculate Laplace transform of the sum H2.61) 

pL P x(s) 


L s (s) — Eexp{—sS 1 } = 


1 - (1 -p)Lpx(s)' 


( 2 . 8 ) 


Let us find transforms of L PX (s) and L s (s ) for a particular case of distribu¬ 
tion of prices Pj, namely, for the case of exponentially distributed prices. In 
this case we have 


L P (s) 


1 

1 + as’ 


where a > 0 is a parameter of exponential distribution. The expression for 
Laplace transform of PjXj is 


r , . y 1 7 tt/- i _7 )s1 _ r (l + l/(ag))r(l + 7 ) 

PX ^ 1 + ask k 11 i T(1 + 7 + l/(as)) 

where T(z) is Euler gamma-function. From (fH2j) we get 

p(r(l + 7 + l/(as)) - T(1 + 7 )r(l + l/(as))) 

Ls(s) = pT(l + 7 + 1 /M) + (1 - p)r(l + 7 )r(l + 1 /(as)) ’ (2-9) 

where the parameters satisfy to the conditions: a > 0 ; 0 < 7 < 1. It is not 
difficult to calculate the limit 

lim (1 - Ls(s))/s 7 = a 7 r(l + 7 )/p. 


In view of this, the distribution of random variable S has heavy (power) tail 
of order 7 . 

Random variable S is not the profit of the short-seller. This is the amount 
he/she received from the sale of assets. Short-seller has to buy the assets to 
give them back. It is necessary to buy J2j=i Xj assets for the price less 
or equal to P*. The profit will be greater or equal to (Pj — P*)Xj. 
Under conditions imposed above, this is a random variable with heavy tailed 
distribution. 
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The list of toy-models may be prolonged, however, it is enough to explain 
multiple ways for tempering such heavy tails distributions. Some of the ways 
will be proposed in the next section. 

3 Tempering procedures 

In this section we propose different variants of tempering procedures. The 
procedures are connected with an attempt to make underlying toy-model 
more precise in a sense. 

3.1 Considerations of Example [2TT 

Example 12.11 has physical sense, and clearly, there are some idealizations in 
the Example. One of idealizations consists in supposition that Brownian 
motion has no drift. The presence of heavy tails in the first passage time 
distribution is connected exactly with this fact. While considering the motion 
with drift, we obtain Inverse Gaussian Distribution as the first passage time 
distribution (see, for example i)- The form of probability density of this 
distribution shows the presence of exponential tail instead of heavy tail of 
Levy distribution. The density of Inverse Gaussian Distribution has form 



(3.1) 


for x > 0. The density of Levy distribution is 



(3.2) 


x > 0. For the case A = a we have 

Pig(x, a, /i) = p L {x, a) exp{a / p}exp{-ax/(2p)}, 

where the term exp{—ax/(2p)} responds for the drift, and exp{<r/p} is a 
normalizing constant. In terms of characteristic functions, this transforma¬ 
tion corresponds to transition from real axis to a line, parallel to it in complex 
plane. This can be seen from their characteristic functions: 













for Inverse Gaussian Distribution (A — a), and 


f L (t) = exp(-v /z 2 Jot) 

for Levy distribution. Transformation from fj J to is exactly “classical” 

tempering procedure (see panam ). Note, that in the considered situation 
the transformation looks absolutely natural, and we have no objection about 
it. 

In a similar way may be considered the case of positive stable distributions 
with index of stability a G (0,1). Namely, let us consider stable random 
variables with Laplace transform 

L st (s) = exp (—As"). 

Multiplying the probability density of this distribution by exp(— ax) and 
making renormalization we obtain new probability density with Laplace 
transform 

Ltst(s) = exp(-A(s + a)") exp(Aa"). 

Corresponding distribution has exponential tail. It is natural to call this 
distribution “tempered” stable distribution. However, this “tempering” has 
no such nice physical interpretation as in the case of a = 1/2, and the 
question on uniqueness and natural character of this procedure is open. 


3.2 Considerations of Example 12.2 

As it has been mentioned above, Example 12.21 is discrete variant of Example 
12.11 Obviously, the distribution with probability generating function 

n*) = (s-s) 

A 

has heavy tail. Because this distribution has the first passage time interpre¬ 
tation, it is natural to propose tempering procedure by introducing a drift 
into random walk. To this aim, it is sufficient to change the distribution by 
the first passage time for the case of non-equal probabilities of moving to the 
right and to the left. For this case, it is known (see m) that the probability 
generating function has the form 


Q(*) 


1 — 1 — 4p(l — p)z 2 

2(1 — p)z 


(3.4) 
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Here p is the probability of moving to the right. We suppose, 1/2 < p < 1. 
In terms of characteristic functions we have 


m 


1 - y/l - e 2it 

e it 


for characteristic function corresponding to (13.3j) . and 


(3.5) 


9(t) 



1 — — exp(2 i(t — ia )) 

exp(*(i — ia) 


(3.6) 


where a = (1/2) log(4p(l — p)), for the case of (j3.4j) . 

We see, that the procedure used in Example 12.II is applied here. In other 
words, classical tempering works nicely in this Example 12.21 as well. 

However, the situation with tempering in Example 12.21 depends on the 
model interpretation. Namely, it may be imagined, that the player (in game 
interpretation of the random walk) has only bounded sum of money to pay 
for each move. The drift is absent. Suppose, the sum of money is M > 0, and 
the player has to pay one unit for each move. Suppose also there is no drift 
in the game (walk). So, we have truncated first passage time distribution. 
Probability generating function now has the form 



(3.7) 

It is obvious, this distribution is not classical tempered variant of that from 
Example 12.21 

Our opinion is that the variant of tempering should be chosen in accor¬ 
dance to the needs of the model correction. The second variant of tempering 
is nor better neither worth than the first one, if there is no information about 
preferable type of correction. Now one can understand, such problem may 
arise also while constructing tempering procedure in Example 12.11 For ex¬ 
ample, we may have time limitation for a particle passes through given level, 
what leads to truncation. 


3.3 Considerations of Example 12.3 

In Example 12.31 we have a product X/l l//2 , where X is zero-mean Gaussian 
random variable, and A is positive stable distribution. Therefore, it is natural 
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to suppose, that tempering procedure for the product has to be based on a 
transformation of stable random variable A. As it was mentioned above, 
we have two variant for the tempering procedure. First one is a classical 
tempering procedure, while the second is a truncation of the distribution of 
A. Third possibility is just a truncation of the distribution of the product 
XA 1 / 2 . Below we show that all three possibilities leads to different results. 

Let us consider these possibilities in details. Suppose that a positive 
random variable A has stable distribution with the index of stability a G 
(0,1) and Laplace transform L(s) = exp{—s Q }. Suppose that a random 
variable X has standard Gaussian distribution. Then the product X ■ A 1 / 2 
has symmetric stable distribution with index of stability 2 a and characteristic 
function 

f{t) = exp{—1 2 “/2"}. 

1. Let start with first tempering procedure of the random variable A. 
Laplace transform of tempered random variable A a has form 

L a (s) = exp{ —(s + a) Q } exp{a a }, 

where a is a parameter of tempering procedure. Characteristic function 
of the product XAd is 

f a (t) = exp{ —(t 2 /2 + a)"} exp{a a }. (3.8) 

The distribution obtained has exponential tails. It converges to initial 
symmetric a-stable distribution as a — > 0. This variant of tempering 
is of an interest because it proposes an application of the procedure to 
one (stable) multiplier, not to the whole product. 

2. The second possibility is to rewrite the product XA l / 2 in the form 
YB where Y is a symmetric stable random variable with the in¬ 
dex /3 G (0, 2) and B is a positive stable random variable with index 
7 G (0,1). In the case when 2a = /hy it is possible to choose the pa¬ 
rameters of Y and B so that XA 1//2 = YB 1 ^. Proposed tempering 
procedure consists now in tempering of B in the way, proposed in item 
03 Obtained product has the tail equivalent to that of the random 
variable Y. Tempering here does not provide exponential tails. Tails 
remains to be heavy, but not so heavy like for initial distribution of 
XA 1 / 2 . 
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3. Third procedure consists in truncation of the random variable A, i.e. 

instead of A we consider Am = min(4, M) for a a fixed constant M > 0. 

1 /2 

Characteristic function of the distribution of XA is 

rM 

E exp-jyCYA^ 2 } = iEexp{—f 2 /24^ 2 } = / exp{— xt 2 /2}dF a (x) : 

Jo 

where F a (x) is distribution function of A. It is clear, that the distri¬ 
bution of Am has Gaussian tails. Corresponding distribution is not 
infinitely divisible. 

Some other variants of tempering are possible as well, but we do not discuss 
them here. 


3.4 Considerations of Example 12.5 

It is more or less clear, that the limit distribution of Z p has heavy tails 
because random number of multipliers u p is large in mean, and u p —> oo in 
probability as n —?• oo. Of course, this situation represents an idealization 
of the reality. It is natural to suppose, that random variable u p is bounded 
with probability 1 by a constant. Therefore, it seems to be natural to change 
geometric random variable v p by bounded random variable /jl Pi m'- 

p(l — n ) k ~ 1 

IP{ft,,M = k} = (3.9) 

for p e (0,1) and integer M > 1. Corresponding probability generating 
function has the form 

pM) (l-(l-p) M )(l-(l-p)z)- 


It is clear, that 


■f* p,M (^) 


pz 


M-> oo 1 — (1 — p)z' 

i.e. the limit distribution in this case is geometric one. However, 

, 2^fc= i* 


p,M (*) 


p —>0 


M 


(3.10) 


(3.11) 


and the limit distribution is uniform over integers from 1 to M. 
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From (I3.10|) and (13. lip is clear, that it is impossible to change the order 
of limits over M and over p. What is a tempered variant of Pareto distribu¬ 
tion? Unfortunately, there is no definite answer to this question. It is possible 
to choose many different possibilities. For example, one can take a sum of 
random number (13.9(1 of arbitrary summands. Because the summands may 
have different distributions, there are infinite many possibilities for the tem¬ 
pering definitions. Another way is the following. One may consider Mellin 
transform of the Pareto distribution as a characteristic function of an expo¬ 
nential distribution. This distribution is geometric stable, and one can make 
classical geometric tempering procedure to get corresponding analogue of ex¬ 
ponential distribution. Now it is necessarily to come back by considering the 
characteristic function as Mellin transform of tempered Pareto distribution. 
Of course, there are many other possibilities to make tempering procedure 
for Example 12.51 


3.5 Considerations of Example 12.6 


It is clear that the reason for the distribution with Laplace transform (12.911 
to have heavy tails is that the distribution of random variables Xj, (j = 
1,2,...) (it is so-called Sibuya distribution) has such tails. The most natural 
way to introduce tempering procedure for this example is to change Sibuya 
distribution by its truncated variant, i.e. instead of probability generating 
function (12.71) use 


Qm(z) 


Ef„(-1)‘ +I 0 


(3.12) 


Of course, corresponding distribution of S has exponential tails. 

As usual, this is not unique way of tempering. We may choose another 
probability generating function 


1 - (1 - azy 
~ 1 - (1 - a)i 


(3.13) 


where a > 0 is tempering parameter. 

Variant ((3. 13(1 is more similar to classical tempering procedure. However, 
the variant (13.121) seems to be more natural one, because short-seller may 
have only a bounded number of the assets. The number of buyers is bounden 
as well. Therefore, they can buy bounded number of assets. 
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3.6 Other ways of tempering 

Classical procedure of tempering is based on exponential truncation of Levy 
measure in Levy-Khinchine representation of stable characteristic function. 
Obviously, it is possible to make some other modifications of this measure. 
For example, it may be changed by a discrete one. Such changes were pro¬ 
posed in [15] in connections with definitions of discrete stable distributions. 
Unfortunately, in [T5j there are no toy-models explaining such way of tem¬ 
pering. 


4 Conclusion 

Morality of the article is that there exist no universal tempering procedure. 
There are many such procedures for a distribution. The choice of tempering 
procedure has to be based on the model (or, at least, on the toy-model) 
leading to this distribution to preserve practical sense of the model, and, in 
the same time, to make more realistic this ideal model. Unfortunately, there 
are almost no toy-models in economics to explain the appearance of heavy 
tailed distributions. It would be nice to have such models to understand 
what type of tempering should be applied. 
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